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The development of sensitive cryogenic light detectors is of primary interest for bolometric experiments 
searching for rare events like dark matter interactions or neutrino-less double beta decay. Thanks to their 
good energy resolution and the natural multiplexed read-out. Kinetic Inductance Detectors (KIDs) are par¬ 
ticularly suitable for this purpose. To efficiently couple KIDs-based light detectors to the large crystals used 
by the most advanced bolometric detectors, active surfaces of several cm 2 are needed. For this reason, we 
are developing phonon-mediated detectors. In this paper we present the results obtained with a prototype 
consisting of four 40 nm thick aluminum resonators patterned on a 2x2 cm 2 silicon chip, and calibrated 
with optical pulses and X-rays. The detector features a noise resolution cte — 154 ± 7eV and an (18±2) % 
efficiency. 


I. INTRODUCTION 

From their first applications in photon detection 1 . Ki¬ 
netic Inductance Detectors (KIDs) became rapidly sub¬ 
ject of several R&D activities in different physics sec¬ 
tors such as astrophysics 2 3 and search of dark matter 
interaction^, proving to be very versatile devices^. In 
superconductors biased with an AC current. Cooper 
pairs oscillate and acquire kinetic inductance. Inter¬ 
actions with energy larger than the binding energy of 
Cooper pairs (2Ao) can break them into quasiparticles, 
producing a variation in the kinetic inductance. The 
superconductor is inserted in a high merit factor LC 
circuit, and the energy of the interactions can be recon¬ 
structed by monitoring the variations in the resonance 
parameters. 

The key features of KIDs consist in the good intrinsic 
energy resolution, and in the possibility of easily tun¬ 
ing each detector to a different resonant frequency. This 
natural aptitude to frequency multiplexed read-out al¬ 
lows operation of a few hundreds of KIDs with two 
cables and one cryogenic amplifier. In addition, since 
the electronics is located at room temperature, the in¬ 
stallation of KIDs arrays would require minor modifi¬ 
cations to the pre-existing cryogenic facilities. Another 
important advantage of these devices is that they are 
operated well below the superconductor critical tem¬ 
perature, where the quasiparticle lifetime and the in¬ 
ternal quality factor saturate, resulting in a negligible 
dependence on the temperature. 
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Light detectors (LDs) based on KIDs could en¬ 
able particle identification, and therefore background 
suppression^ 7 in bolometric experiments searching for 
neutrino-less double beta decay (Ov/3/3) or dark matter 
interactions, such as CUOREP" and LUCIFEKp. 

The present limit of KIDs in this kind of applica¬ 
tion resides in their small active surface (few mm 2 ). 
Macro-bolometers used by CUORE and LUCIFER fea¬ 
ture a surface of several cm 2 , thus an effective read¬ 
out of the light emitted by these crystals demands for 
LDs with large active area. Obtaining a similar surface 
with KIDs would require hundreds of pixels for each 
LD, which is not a realistic option for experiments with 
~ 1000 detectors. This problem can be overcome fol¬ 
lowing the phonon-mediated approach developed for 
Neutron Tansmutation Doped Ge thermistors 10 and for 
Transition Edge Sensors 11 , and recently proposed also 
for KIDsPEE]. 

The CALDER project 14 aims to demonstrate that such 
an approach allows to develop large-area LD with noise 
energy resolution <20 eV RMS, multiplexed read-out, 
high reproducibility and scalability. 

In this paper we present the results of a CALDER 
development prototype, implementing four aluminum 
resonators. We perform an absolute calibration of the 
energy absorbed by the KIDs and derive the efficiency 
of the detector. We characterize the response with pho¬ 
tons at 400 nm and X-rays from a 57 Co source. Finally, 
we analyze and evaluate the energy resolution of the 
detector. 
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II. DETECTOR FABRICATION 


DATA ANALYSIS 


The detectors are fabricated at CNR IFN on high 
quality, 300/im thick, high resistivity (>10 kfixcm) 
Si(100) substrates. The four lumped-element resonators 
(in the following KID-1, KID-2, KID-3 and KID-4, or¬ 
dered according to the position on the chip) are pat¬ 
terned by electron beam lithography on a single 40 nm 
thick A1 film deposited using electron-gun evapora¬ 
tor. The active area of the single pixel consists of 
an inductive meander made of 14 connected strips of 
80 }tm x 2 mm. The meander is closed with a capacitor 
made of 5 interdigitated fingers of 1.2 mm x 50 pm, to 
ensure the uniformity of the current across the induc¬ 
tor. The resonant frequency of each resonator is varied 
by cutting the last finger of the capacitor. As shown in 
Figure [lj the chip is assembled in a copper structure 
using PTFE supports with total contact area of about 
3 mm 2 . The other side of the holder (not shown) is 
covered with a copper collimator hosting a 57 Co cali¬ 
bration source (peaks at 6.4 and 14.4 keV) and an opti¬ 
cal fiber coupled to a room-temperature LED, that pro¬ 
duces pulses at 400 nm. The source and the fiber are 
placed on the back of the substrate to avoid direct il¬ 
lumination of the resonators. The optical fiber points 
to the center of the substrate, while the 57 Co source is 
located nearby KID-4. 



FIG. 1. The four A1 KIDs deposited on the 2x2cm 2 Si sub¬ 
strate. The chip is assembled in a copper structure using 
PTFE supports and illuminated from the back with a colli¬ 
mated 57 Co source and an optical fiber. 

The copper holder is thermally anchored to the cold¬ 
est point of a 3 FIe/ 4 FIe dilution refrigerator with base 
temperature of 10 mK. 

The output signal is fed into a CITLF4 SiGe cryogenic 
low noise amplifier 15 operated at 4K. A detailed de¬ 
scription of the chip design, the cryogenic setup of our 
laboratory at Sapienza University in Rome, the room- 
temperature electronics and the acquisition software 
can be found in references 141617 . 


A typical data collection consists in acquiring the 
complex transmission (S21) for a frequency sweep 
around the resonances (see Figure 2j, and fitting the 
resonance circles in order to extract the quality factor Q, 
the coupling quality factor Q c and the internal quality 
factor Q, using the method described in references 1819 . 



FIG. 2. Amplitude of the line transmission (S21) for a VNA 
frequency sweep around the resonances. Inset: fit of the reso¬ 
nance circle of KID-3 in the working point; the green marker 
indicates the resonant frequency. 


Resonators are designed to be over-coupled, thus the 
total quality factors Q (reported in Table fl| are entirely 
dominated by the coupling quality factors Q c . Q,- 
values differ from the design value of 8000, likely due 
to the presence of parasitic slotline modes or standing 
waves between the detector and the cryogenic amplifier, 
or inter-resonator coupling resulting in coupled oscilla¬ 
tion mode^. Q; is above 150 x 10 3 , but the values of Q, 
limit the accuracy of the estimation. 


TABLE I. Resonant frequency /o, quality factor Q, optimal 
(off-resonance) microwave power at the amplifier input P,„ 
and quasiparticles recombination time T pp . Errors on fg and 
P,„ are negligible, while errors on Q and T^p are dominated 
by fit systematics and are lower than 10%. 



/o 

Q 

P 

1 in 

T q V 


[GHz] 

[xlO 3 ] 

[dBm] 

[?«] 

KID-1 

2.675 

6 

-63 

218 

KID-2 

2.689 

18 

-64 

228 

KID-3 

2.731 

8 

-66 

243 

KID-4 

2.746 

35 

-72 

233 


When the trigger of any of the resonators fires, we 
acquire a 2 ms long time window for the real and imag¬ 
inary parts of S21 for each resonator (I and Q) with a 
sampling frequency of 500kFIz. I and Q variations are 
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then converted into changes in phase (Sep) and ampli¬ 
tude relative to the center of the resonance loop (blue 
marker in the circle reported in Figure |2|. In the fol¬ 
lowing analysis we use only the Sip signal, as for this 
detector it is from 6 to 10 times larger than the ampli¬ 
tude one, depending on the KID. 

To determine the optimal microwave power, we eval¬ 
uate the signal-to-noise ratio scanning from -80 dBm to 
-50 dBm. Increasing the input power produces a reduc¬ 
tion of the noise contribution from the amplifier but, on 
the other hand, decreases the quasiparticles recombina¬ 
tion time Tqp and, as a consequence, the signal integra¬ 
tion length. The microwave power that optimizes the 
signal to noise ratio for each resonator is reported in 
Table 0 

The average noise power spectrum is reported in Fig¬ 
ure [3] for phase (continuous line) and amplitude (dot¬ 
ted line) read-out of each resonator. The flat noise ob¬ 
served in the amplitude read-out and in the high fre¬ 
quency region of the phase read-out, is consistent with 
the noise temperature of the amplifier (Tjv ~ 7 K). The 
low-frequency region of the phase spectra is dominated 
by another noise source, whose origin is not clear yet. 
It is not ascribable to two-level system noise, as it does 
not depend on temperature or microwave power. Fur¬ 
thermore, the presence of a mu-metal shield around 
the cryostat should guarantee an efficient suppression 
of noise due to static or low-frequency magnetic fields. 
Since a fraction of this noise is found to be correlated, 
it could be caused by frequency jitters in the read-out. 



FIG. 3. Average noise power spectrum for phase (continu¬ 
ous line) and amplitude (dotted line) read-out. On top of the 
white noise from the amplifier, the phase noise exhibits an 
extra contribution at low frequency. 



FIG. 4. Response of the four resonators to 15 keV pulses pro¬ 
duced by the optical fiber, placed in the proximity of KID- 
2 and KID-3. The responses are obtained averaging many 
pulses to reduce the random noise. 


random noise contributions. 

We fit the pulses with a model that includes the time 
constant of the low-pass filter, the ring-time of the res¬ 
onator (r r = Q / (nfo )) and two free parameters: a rise¬ 
time, which is related to the spread in the arrival time 
of phonons, and a decay time. As expected, the rise¬ 
time depends on the distance between the resonator 
and the optical fiber, and ranges from 2 fis (KID-2 and 
KID-3) to 10 [is (KID-1) and 17 i/s (KID-4). The decay 
time becomes faster increasing the microwave power or 
the temperature, and for this reason it is identified as 
Tqp. This time constant does not depend on the energy 
of the optical pulses in the scan range (0.7-25 keV), and 
its value is reported in Table [I] for each resonator. 


IV. ENERGY CALIBRATION AND EFFICIENCY 

The energy E absorbed in a resonator creates a num¬ 
ber of quasiparticles SNqp — i]E/ Aq, where // is the de¬ 
tection efficiency. The variation SNqp produces a linear 
shift of the resonant frequency / from the equilibrium 
one /(P: 


E = £±SN, p = 


A o 

V 



f~fo \ 
fo ) 


( 1 ) 


The high frequency noise in the acquired waveforms 
is rejected off-line using a software low-pass filter. In 
order to avoid distortions in the rise-time of the pulses, 
the cut-off frequency is set at 100 kHz (t cu t- 0 ff ~ 
1.6 ps). Finally, the waveforms are processed with the 
optimum filter 21 22 , which includes a resonator-specific 
rise time (see below). The results are not highly sensi¬ 
tive to the choice of rise time. 

In Figure]?] we show the typical response of the four 
resonators to the interaction of 15keV optical pulses, 
obtained by averaging several pulses to suppress the 


where po — a.S 2 (f, T)/ANoVAq. The parameter a 
is the fraction of the total inductance due to kinetic 
inductance. No is the single spin density of states 
(1.72x 10 10 eV -1 pm -3 ) and S 2 {f,T) is a slow function 
of the temperature and of the resonant frequency that 
relates the phase variation of the complex conductiv¬ 
ity to Cooper pairs breaking. In our working condi¬ 
tions, S 2 {f,T) is measured to be 2.3-2.6 depending on 
the resonator. The active volume of the resonator V 
is calculated by correcting the volume of the inductor 
(96500 pm 3 ) for the average variation of current density 
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(90%) evaluated from SONNET simulation^®. 

We calculate Ao and a by measuring the variations 
of the resonant frequency (/ — /o)//o as a function of 
the temperature, and fitting these data to the Mattis- 
Bardeen theory according to the procedure described 
by Popel 24 and adopted by Gao et al , 25 . The relation¬ 
ship between (/ — /o)//o and T, indeed, depends only 
on Aq and a, that are found to be the same for all the 
resonators. The average values are Ao = (201 ± 6) peV 
and a. — (5.8 ± 0.5)%. 

We make a second, independent measurement of Ao 
by illuminating the chip with a Millimeter-Wave source. 
We monitor phase, amplitude and resonant frequency 
of the resonators while increasing the frequency of the 
source from 75 to 110 GHz. The minimum source fre¬ 
quency that produces significant variations of the res¬ 
onators features is v m — 95.5 GHz. This value corre¬ 
sponds to Ao = hv m /2 = (197 ± 5) ;<eV, in full agree¬ 
ment with the previous value. With these parameters, 
po turns out to be (1.2 ± 0.1) x 10~ 13 for all the res¬ 
onators. 

Finally, equation [l] can be modified in order to con¬ 
vert the frequency response (/ — /o)//o into the mea¬ 
sured phase variation 8(p recalling that, for Sep mea¬ 
sured in the center of the resonance loop, (/ — fo)/fo — 
Scp/AQ: 


E = 


A o 


4Q>/Po 


5<p. 


( 2 ) 


To estimate 6<p we could use directly the amplitude of 
the phase signal. However, to account for the observed 
small differences in the pulse shapes induced by the ar¬ 
rival time of phonons, we evaluate Sep by integrating the 
time development of the pulse and dividing the result 
by r , qp . 

The energy-calibration function in Eq. [2] is applied 
to the single resonators to obtain the energy spectra of 
the 57 Co source (reported in Figure [ 5 ] (left)) and of the 
optical pulses. To infer the detection efficiency, we fit 
the position of the high energy peak produced by the 
57 Co source, and we divide the result by the nominal 
energy of the line (see Table pi. Applying this proce¬ 
dure to the lower energy peak (6.4 keV) gives consistent 
results. The estimation of the efficiency on X-rays al¬ 
lows us to neglect uncertainties in the energy deposited, 
which could be significant when dealing with optical 
pulses. To evaluate the total efficiency of the detec¬ 
tor, we sum event-by-event the calibrated energies of 
each KID (see Figure [5] (right)) and fit the position of 
the 14.4 keV peak, obtaining i]sum — 2.6 keV/14.4keV = 
(18±2) %. The error on the efficiency is dominated by 
the error on Tqp introduced by the absolute calibration 

(Eq# 

We apply to optical pulses the analysis technique 
validated on X-rays. The energy of the optical pulses 
was previously calibrated with a PMT at room temper¬ 
ature, and corrected for the PMT quantum efficiency. 


xlO 3 xlO 3 




0 0.2 0.4 0.6 0.8 1 1.2 

Absorbed Energy [keV] 


1 2 3 4 5 6 7 8 

Summed Absorbed Energy [keV] 


FIG. 5. Left: low energy spectrum of the energy absorbed in 
KID-1 (red, continuous line), KID-2 (gray, filled) KID-3 (blue, 
dotted line) and KID-4 (green, dashed line) calibrated with 
the function in Eq. |2j each resonator absorbs a different en¬ 
ergy fraction of the X-rays from the 57 Co source depending on 
its position. Right: event by event sum of the the calibrated 
energies reported in the left plot; the fit of the 14.4 keV peak 
produced by the 57 Co source is shown. The energy spectra 
are not corrected for the detection efficiency. 


for the reflectivity of silicon 26 , and for the geometrical 
efficiency, evaluated through a Monte Carlo simulation 
based on the Litrani software!- 27 . We observe that the 
detector response is linear in the scan range. The ener¬ 
gies of the optical pulses evaluated using the absolute 
calibration (Eq. [2j and corrected for TjsuMr differ by less 
than 5% from what computed with the simulation. 

In Table III] we report the noise resolution erg, evalu¬ 
ated at the detector baseline, and the intrinsic noise of 
each resonator, calculated by scaling the resolution oy 
for the single KID efficiency tj. Summing the response 
of the four KIDs on an event by event basis, we obtain 
a global energy resolution o~e,SUM of 154±7eV, which 
corresponds to an intrinsic resolution of cte,SUM X ?/ = 
154 eV x 0.18 = 27.7eV. This resolution is worse than 
the one expected if we were dominated by the amplifier 
noise that, for each resonator, can be computed as 14 : 


amp _ AqNoV Mq 2 

E “ rjoLQS 2 (J,T)^ U 

The discrepancy between the measured values c7j and 
Up'' 11 ’ can be interpreted as due to the excess low 
frequency noise in the phase readout. This is sup¬ 
ported by the fact that (Je matches closely c r f SD , which 
is the integral of the optimally-filtered noise power 
spectral density. Without this noise, we would have 
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TABLE II. Measured noise resolution (cr E ), expected noise res¬ 
olution in the limit in which the amplifier noise dominates 
(Cg and integral of the optimally-filtered noise power spec¬ 
tral density (a^ SD ); pixel efficiency ;/ and intrinsic resolution 
of the resonator (cr E x ;/). In the last line we report the global 
performance of the detector. 



O-E 

^amp 

PSD 

V 

<t e X 1] 


[eV] 

[eV] 

[eV] 

[%] 

[eV] 

KID-1 

400 

147 

484 

3.1±0.4 

12.4 

KID-2 

262 

60 

253 

3.4±0.4 

8.9 

KID-3 

205 

75 

233 

6.1±0.7 

12.5 

KID-4 

204 

50 

184 

5.5±0.6 

11.2 

SUM 

154 

42 


18±2 

27.7 


expected a global energy resolution of a ™'sum ~ 
y/U47 X Vi) z )/’1SUM = 42eV. 


V. CONCLUSIONS 

In this paper we presented the results obtained with 
a 4-pixels A1 array on Si substrate. The comparison of 
the absolute calibration of the detector with the energy 
produced by a calibrated LED+optical system proved 
the reliability in reconstructing the energy of optical 
pulses from 0 to 25keV. We derived an efficiency of 
(18±2) %, that can be improved by increasing the active 
surface of the resonators. The overall baseline resolu¬ 
tion of 154±7 eV is already competitive with some of 
the commonly used cryogenic light detectors, and can 
be further improved by suppressing the noise sources 
and, eventually, by using more sensitive superconduc¬ 
tors. 
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